Method for determining the position of the position indicator of a position measuring system

ABSTRACT

The invention concerns a method for determining a position φ of a position indicator of a position measuring system, in which the sensor signals a 1  to a r  of a number of r sensors are recorded as a signal vector {right arrow over (a)}=(a 1 , . . . , a r ), a measurement vector {right arrow over (q)} is formed according to {right arrow over (q)}=g({right arrow over (a)}), a 2-component vector {right arrow over (p)} is calculated according to {right arrow over (p)}=M·{right arrow over (q)}, wherein M is a 2×n matrix, and the position φ is determined by means of a predetermined function ƒ({right arrow over (p)}) to φ=ƒ({right arrow over (p)}), wherein the function ƒ({right arrow over (p)})=ƒ(p 1 ,p 2 ) is based on the equations p 1 =cos(φ) and p 2 =sin(φ). The matrix M is determined in a calibration phase in such a manner, that the transformation {right arrow over (p)}=M·{right arrow over (q)} is a linear mapping, which maps a number m of measurement vectors {right arrow over (q)} k  on 2-component vectors {right arrow over (p)} k , wherein the tips of the vectors {right arrow over (p)} k  substantially lie on a circle or on a circular arc. The method allows the elimination of external fields.

PRIORITY CLAIM

Applicant hereby claims foreign priority under 35 U.S.C § 119 from Swiss Patent Application No. 01207/16 filed on Sep. 16, 2016 and from European Patent Application No. 17167135.7 filed on Apr. 19, 2017, the disclosures of which are herein incorporated by reference.

FIELD OF THE INVENTION

The invention concerns methods for determining the position of the position indicator of a position measuring system, in particular a method for determining the position of a position indicator of a linear position measuring system or of the position indicator of an angle measurement system.

BACKGROUND OF THE INVENTION

Various sensor systems for measuring the position φ of a rotating or linearly moving position indicator are known, in which the information to be obtained is represented by sine and cosine shaped varying raw signals p₁=cos(φ) and p₂=sin(φ). Possible implementations include e.g. the measurement of components of the induction field of a permanent magnet at suitable locations in space. The searched position is then determined from p₁, p₂ e.g. via the function φ=a tan 2(p₂,p₁) by means of the function a tan 2 known from the C-programming language. These prerequisites severely restrict possible arrangements of sensor elements and position indicator, and in particular require a high degree of accuracy in their production or the relative positioning of the sensor chip and the position indicator, respectively.

A fundamental difficulty with these methods is that a non-ideal arrangement of the position indicator and the sensor chip leads to distortions of the raw signals, in particular that these no longer have the pure sine and cosine shape and thus lead to an inaccurate determination of the position φ. There exist methods to improve the raw signal quality, e.g. by using offset and amplitude corrections. A method for a sensor arrangement with two sensors known from WO 2005124286 is based on a description of the measuring points as ellipses and requires specific calibration measurements, which provide parameters for the corrections and which can be realized in practice only with considerable effort. This method requires great care in the design of the sensor arrangement and especially in the adjustment of the permanent magnet, so that the measured values lie sufficiently precisely on an ellipse. If the actual rotational axis does not pass through the magnetically defined center of the permanent magnet, the correction can even cause the opposite effect of an increase in the angular error.

In many instances, a position sensing system is responsive to external fields (e.g., the earth's magnetic field in the case of a position indicator based on permanent magnets). Without countermeasures, these additional field components cannot be distinguished from the useful field so that the external fields restrict the accuracy of the position measuring system.

A fundamental challenge for position measuring systems is their long-term stability, i.e., an initially sufficiently accurate calibration changes over time, e.g. by mechanical displacement of the components relative to one another or by drifting of electronic sensor properties. However, the measurement quality of systems during normal operation is not ascertainable according to the state of the art, in particular this cannot take place intrinsically, i.e. without reference measurements from outside. However, an intrinsic method would be of great importance for fault-tolerant systems. In addition, a recalibration of a position measuring system during normal operation would be very valuable in order to decisively improve the long-term stability.

SUMMARY OF THE INVENTION

The method according to the invention for determining the position of a position indicator is based on the finding that the sensor measurement values supplied by the sensors can be represented as a measurement vector, i.e. as a point in a vector space, the dimension of which corresponds to the number of sensors, and that the measurement vectors generated during the movement of the position indicator (rotation about an axis of rotation or displacement along an axis) lie on a path called orbit in this vector space.

This orbit can—under ideal conditions—be projected onto a circle or a circular arc by means of a linear mapping, i.e. each position of the position indicator corresponds to a point on the orbit in the abstract space of the signal vectors which is unambiguously assigned to a point on the circle or arc. Non-ideal conditions, for example non-ideal properties of the position indicator, which can be particularly the case with magnets, or an insufficient number of sensors or the use of non-optimally selected calibration positions, cause deviations from the circular shape. The linear mapping can be represented as a matrix-vector operation and can be determined using standard methods of linear algebra. Each measurement vector is projected onto a 2-component vector, the tip of which approximates a circular arc, wherein the direction of the 2-component vector corresponds to the position of the position indicator. The linear mapping is determined in a calibration phase, which includes the detection of measurement vectors and the decomposition of a matrix formed by the measurement vectors into two or typically three matrices. The representation of the measured values provided by the sensors as a vector also makes it possible to reduce the influence of external fields, in particular magnetic external fields, by means of a further linear mapping.

The method according to the invention can be used universally. The arrangement of the sensors is virtually arbitrary. The sensors may even detect different components of the field generated by the position indicator. The determination of the matrix can be carried out using measurements which only have to fulfill an absolute minimum of prerequisites. The method is particularly effective without a plurality of angular reference measurements, except the determination of the zero position, which however requires only a single reference measurement. These minimum requirements for the measurements needed for the calibration allow inter alia that a recalibration of the position measuring system is possible during normal operation, or that a factory-based, initial calibration for a given position measuring system can be iteratively improved by means of a bootstrapping method.

Further, in implementations which provide redundant information to the position of the position indicator (e.g., because more than the minimum number of single sensors is installed), one or more of the following additional tasks can be solved:

-   -   error detection, i.e. intrinsically detecting a drift and         triggering a predetermined action, for example, a corresponding         warning or recalibration.     -   external field compensation by suppressing external field         influences on the measuring signals.     -   elimination of errors caused by an offset between the actual         axis of rotation of the position indicator and the ideal axis of         rotation.     -   elimination of errors caused by irregularities of the magnet         when the position indicator is a magnet.     -   improvement of the signal-to-noise ratio if more than two         sensors are used, since then a plurality of individual sensor         values enter into the determination of the position and thus an         averaging takes place. In particular, difficulties during the         zero crossing of the sine or cosine signal of conventional         systems can be avoided.

Furthermore, the method for determining the position of a position indicator according to the invention can be designed such that a number of r=2+b sensors are used, and that a projection matrix is used to project the influence of a number of b different external field configurations out of the signal vectors, wherein the projection matrix simultaneously produces 2-component vectors whose tips lie essentially on a circle or circular arc.

BRIEF DESCRIPTION OF THE DRAWING FIGURES

The accompanying drawings, which are incorporated into and constitute a part of this specification, illustrate one or more embodiments of the present invention and, together with the detailed description, serve to explain the principles and implementations of the invention. The figures are drawn schematically and not to scale. In the drawings:

FIG. 1 shows an angle measurement system,

FIG. 2 shows a linear position measuring system,

FIG. 3 shows a resolver, and

FIGS. 4 to 10 illustrate various aspects of the method according to the invention.

DETAILED DESCRIPTION OF THE INVENTION

FIG. 1 shows a position measuring system which has a magnet 2 rotatable about a predetermined axis of rotation 1 and a sensor chip 3. The rotational position of the magnet 2 is the position φ to be determined by means of the sensor chip 3. The sensor chip 3 is shown enlarged on the right-hand side of FIG. 1. The sensor chip 3 comprises a predetermined number of 1 to r sensors 4 and electronic circuits for the operation of the sensors 4, etc. The sensors 4 are, for example, Hall sensors. Such a sensor 4 can also be a cluster of several individual sensors which provide a common output signal. The sensor chip 3 can also contain magnetic field concentrators, which amplify the magnetic field generated by the magnet 2 at the locations of the sensors 4. The number r of the sensors 4 is at least r=2. The sensors 4 each output a sensor signal. The sensor signals are designated herein as a₁ to a_(r). The index j of the sensor signal a_(j) therefore designates the number of the respective sensor 4.

FIG. 2 shows a position measuring system constituted by a magnet 2 and a sensor chip 3 which are displaceable relative to one another along a linear axis x. The sensor chip 3 comprises a predetermined number of 1 to r sensors 4 and electronic circuits for the operation of the sensors 4, etc. The sensors 4 are, for example, Hall sensors. The magnet 2 consists of alternately oppositely magnetized regions which produce a periodically changing magnetic field at the locations of the sensors 4 during a displacement along the axis x. Die x-position of the magnet 2 is the position φ to be determined by means of the sensor chip 3. It is x=D/2π·φ, wherein D is a distance defined by the magnet 2.

The sensors 4 of these position measuring systems provide, as a result of the periodically changing magnetic field, harmonic output signals. The position measuring system of FIG. 1 is an angle measuring system with the rotatable magnet 2 as a position indicator. The position measuring system of FIG. 2 is a linear measuring system with the displaceable magnet 2 as a position indicator. In both position measuring systems, the range of values of the position φ extends over a range from 0 to 2π, which corresponds either to the angular range from 0 to 360° or to a linear range of 0 to D, wherein the distance D results from the fact, that for all j=1, . . . , r sensor signals a_(j) is a_(j)(x+D)=a_(j)(x).

FIG. 3 shows a position measuring system, which is formed from a rotor coil R and from the stator coils S_(i), i=1, . . . , r with r≧2. The rotor coil R is rotatable about an axis of rotation and is fed with a sinusoidal current which generates a magnetic field. The magnetic field induces in the stator coils S_(i), i=1, . . . , r a voltage whose amplitude is modulated as a function of the rotational angle r of the rotor coil R. By means of demodulation, sensor signals a_(i)(t), i=1 . . . , r, which vary harmonically with the angle of rotation φ, are generated from the stator coil signals, from which the instantaneous rotational angular position φ(t) of the rotor coil R can be calculated. The signal vector a_(i)(t) consists of the envelope of the signal of the stator coil S_(i) for all i=1, . . . , r.

Position measuring systems can also be based on the inductive principle, in which at least one transmitting coil and at least one receiving coil are spatially fixed, and the position indicator is a geometrically structured, electrically conductive disk which is rotatable about an axis. It is also possible to use the same coil for the transmission and the reception process. The coils can be arranged and/or connected to one another on several levels. The position measuring system can have two or more coils. The use of several coils is particularly helpful to detect the full 360° rotation. The position indicator changes the inductance of the transmitting coil and the receiving coil(s). The transmitting coils are fed with alternating currents, preferably sinusoidal currents, whereby the magnetic fields generate induction voltages in the receiving coils which because of the position-dependent eddy currents generated in the disk have harmonically varying amplitudes. By means of demodulation, the signals for the method according to the invention for determining the position of the position indicator are obtained.

The method according to the invention can also be used in optical position measuring systems in which two or more optical sensors (e.g. phototransistors, photodiodes, etc.) detect in dependence of the position of an optically active position indicator the brightness, polarization, phase shift, time of flight or other optical property variations. For this purpose, the position indicator reflects or modulates a light beam and thus causes the harmonic signals used by the method. Instead of the optical field, such a position meter can also be equipped with radar or terahertz radiation sensors, or corresponding position indicators. In addition, sound waves, in particular ultrasound can also be used instead of optical and electromagnetic waves, and a corresponding position measuring system can be constructed.

The method according to the invention can also be used in electrostatic position measuring systems in which a rotatably mounted geometrically structured electrode as position indicator together with stator electrodes forms two or more capacitors with position-dependent capacitances. Capacitive signals that vary with the position are generated and detected. In an electrostatic system, the position indicator can also be realized by a structured dielectric material which modulates the capacitance of two or more spatially fixed electrode pairs harmonically with the position.

The method according to the invention also makes it possible to separate the useful signal or the useful signals from interfering signals. An example is a magnetic rotary encoder on an electric motor. In this case, the magnetic field of the magnetic angle sensor can be separated cleanly from the rotating fields of the electric motor. Thus, the position can be determined without interference. This also applies to optical, electro-magnetic, electrical, magnetic, and acoustic position determinations.

All these physical principles can be used both for the construction of a position measuring system with a rotatably mounted position indicator as well as with a linearly displaceably mounted position indicator.

The invention relates to a method by which the position φ can be calculated from the measured values delivered by the sensors 4 in virtually any position of the sensor chip 3 relative to the axis of rotation 1 or to the linear axis x. The method may be a simple method, e.g. a method without offset correction and without external field correction, or a method with offset correction or a method with external field correction or a method with offset and external field correction. The method comprises, in all cases, a calibration phase and the normal operation.

The offset correction takes place by means of an offset vector {right arrow over (a)}₀. The external field correction takes place by means of a projection matrix P. The offset and external field correction takes place by means of the offset vector {right arrow over (a)}₀ and the projection matrix P.

The offset vector {right arrow over (a)}₀ and the projection matrix P are determined, if necessary, in the calibration phase. In the calibration phase, a matrix M is determined in each case.

Normal Operation

In normal operation, the position φ of the position indicator is calculated by the steps:

-   -   A) Detecting the sensor signals a₁ to a_(r) of the sensors.         -   The sensor signals a₁ to a_(r) form an r-tuple (a₁, . . . ,             a_(r)), which is interpreted as a signal vector {right arrow             over (a)}=(a₁, . . . , a_(r)) in an r-dimensional signal             vector space.     -   B) Forming a measurement vector {right arrow over (q)} according         to a predetermined linear function to {right arrow over         (q)}=g({right arrow over (a)}).     -   C) Calculating a vector {right arrow over (p)}=M·{right arrow         over (q)}.         -   {right arrow over (p)}=(p₁,p₂) is a 2-component vector with             the components p₁ and p₂ and M is a 2×n matrix.

${\overset{\rightarrow}{p} = {{{M \cdot \overset{\rightarrow}{q}}\mspace{14mu} {can}\mspace{14mu} {be}\mspace{14mu} {written}\mspace{14mu} {as}\mspace{14mu} \begin{pmatrix} p_{1} \\ p_{2} \end{pmatrix}} = {\begin{pmatrix} M_{11} & \ldots & M_{1\; n} \\ M_{21} & \ldots & M_{2\; n} \end{pmatrix} \cdot \begin{pmatrix} q_{1} \\ \vdots \\ q_{n} \end{pmatrix}}}},$

-   -   D) Determining the position φ, from the components p₁ and p₂ of         the vector {right arrow over (p)} by means of a predetermined         function ƒ({right arrow over (p)}) to φ=ƒ({right arrow over         (p)})=ƒ(p₁,p₂).

Preferred examples for the linear function g({right arrow over (a)}) are:

-   -   1) g({right arrow over (a)})={right arrow over (a)}, i.e. {right         arrow over (q)}={right arrow over (a)}, if neither an offset         correction nor an external field correction is performed,     -   2) g({right arrow over (a)})={right arrow over (a)}−{right arrow         over (a)}₀, i.e. {right arrow over (q)}={right arrow over         (a)}−{right arrow over (a)}₀, if (only) an offset correction is         performed,     -   3) g({right arrow over (a)})=P·{right arrow over (a)}, i.e.         {right arrow over (q)}=P·{right arrow over (a)}, if (only) an         external field correction is performed, or     -   4) g({right arrow over (a)})=P·({right arrow over (a)}−{right         arrow over (a)}₀), i.e. {right arrow over (q)}=P·({right arrow         over (a)}−{right arrow over (a)}₀), if an offset and external         field correction is performed.         wherein {right arrow over (a)}₀ in designates an offset vector         and P a projection matrix. The measurement vectors {right arrow         over (q)} therefore have a number, n, of components, with 2≦n≦r.

FIG. 4 illustrates the relationship of positions φ_(k) of the position indicator and the vectors {right arrow over (a)}_(k), which are formed from the signals of the r sensors, for k=1, 2 to m. The positions φ_(k), which lie in the range from 0 to 2π, are shown here as angles, as it is the case with an angle measuring system. However, the positions φ_(k) of the position indicator can also lie on a linear axis. The arrows illustrate the positions φ_(k) of the position indicator, a signal vector {right arrow over (a)}_(k) is assigned to each position φ_(k).

FIG. 5 illustrates for the example n=3 the linear mapping {right arrow over (p)}=M·{right arrow over (q)}, which maps the measurement vectors {right arrow over (q)}, which lie on an orbit 5 in a n-dimensional vector space, by means of the matrix M on 2-component vectors {right arrow over (p)}, the tips of which lie on a circle 6. The left side illustrates the space of the measurement vectors {right arrow over (q)}, the right side shows the plane with the circle 6 and the 2-component vectors {right arrow over (p)}. By way of example, the mapping from {right arrow over (q)}₃ on {right arrow over (p)}₃ is indicated by an arrow. A vector {right arrow over (p)}_(k) is assigned to each measurement vector {right arrow over (q)}_(k). The direction of the vector {right arrow over (p)}_(k) includes an angle with a reference direction which corresponds to the position φ_(k) of the position indicator.

Temperature Compensation

The sensor signals may be temperature-dependent. The temperature dependence can be compensated in different ways, for example, in that the temperature is measured and the sensor signals as well as the sensor offsets are temperature-compensated. Alternatively, the above process can be carried out at various temperatures and the relevant quantities, in particular the matrix M, can be a temperature-dependent matrix. The matrix M is determined in this case, for example, for a predetermined number of temperature support points and the temperature compensation is effected as follows:

-   -   1. Detecting the temperature and interpolating the matrix M for         the measured temperature from the matrices M at the support         points.     -   2. Applying the interpolated matrix M to the actual measurement         vectors.

If an offset compensation and/or an external field compensation is carried out, also the offset vector {right arrow over (a)}₀ and given the case the projection matrix P are determined in the same way by interpolation of the corresponding objects at the temperature support points.

This variant offers the advantage that in this way both temperature effects which influence the mechanical position of the position indicator relative to the sensor chip as well as temperature-dependent effects in the sensors can be compensated for.

The vector {right arrow over (p)} contains the complete information about the position φ to be determined. For many applications, the relationship between the vector {right arrow over (p)} and the position φ is given by the equations p₁=cos(φ) and p₂=sin(φ). The predetermined function ƒ({right arrow over (p)})=ƒ(p₁,p₂) is a function, which is based on the equations p₁=cos(φ) and p₂=sin(φ). It is, in particular, a function, which determines the position φ on the basis of the equation φ=aucustangens(p₂/p₁) and the signs of the values p₁ and p₂. The function ƒ({right arrow over (p)}) is for example the function ƒ({right arrow over (p)})=a tan 2(p₂,p₁), wherein a tan 2(p₂,p₁) is a function known from many programming languages which provides the correct angle in the range of 0 to 2π.

Calibration Phase

Since the offset vector influences the determination of the projection matrix P as well as the determination of the matrix M, and the projection matrix P influences the determination of the matrix M, the following sequence must be observed in the calibration phase:

-   -   1. Determination of the offset vector {right arrow over (a)}₀,         if an offset correction is desired,     -   2. Determination of the projection matrix P, if a correction of         external fields is desired,     -   3. Determination of the matrix M.         These steps are subsequently explained in detail.

Determination of the Offset Vector

The offset vector {right arrow over (a)}₀ can be determined in a conventional, known manner. In doing so, the offset of each of the sensors a_(0j) with j=1 to r is determined and the entirety of the offsets is stored as offset vector {right arrow over (a)}₀=(a₀₁, . . . , a_(0r)).

The offset vector {right arrow over (a)}₀ can also be determined from the signal vectors, if the positions belonging to the signal vectors are equidistributed. This occurs by:

-   -   1. Forming the mean values ā₁ to ā_(r) of the single sensor         signals according to ā_(l)=1/mΣ_(k=1) ^(m)a_(kl).     -   2. Forming the offset vector {right arrow over (a)}₀=(ā₁, . . .         , ā_(r)).

Determination of the Projection Matrix P

The signal vector {right arrow over (a)}=(a₁, . . . , a_(r)), which form the orbit, are generally located in a subspace of the entire r-dimensional vector space of the r-tuples. Influences of external fields can be eliminated in this case by the application of a projection matrix P. The determination of the projection matrix P occurs, preferably in absence of the position indicator, with the following steps. However, it can also be performed in presence of the position indicator.

-   -   1. In absence of the position indicator step 1 comprises the         substeps 1.1 to 1.2, which are executed for a predetermined         number of different external fields v=1 to b≧1, wherein the         number b fulfills the condition b≦r−2.         -   1.1 Applying the external field.             -   The external field is for example a homogeneous magnetic                 field, which points in an arbitrary space direction. The                 external field can also be an inhomogeneous external                 field, to which the position measuring system is                 exposed. In order to eliminate the influence of external                 fields, additional sensors 4 are required. The minimum                 number of sensors 4, which is necessary to eliminate a                 number of b different external fields, is according to                 the above equation r=b+2.         -   1.2 Detecting the sensor signals a₁ to a_(r) and forming an             external field vector {right arrow over (ƒ)}_(v)={right             arrow over (a)}=(a₁, . . . , a_(r))             -   or given the case of an offset compensated external                 field vector, respectively, {right arrow over                 (ƒ)}_(v)={right arrow over (a)}−{right arrow over (a)}₀.                 or     -   1. In presence of the position indicator step 1 comprises the         substeps 1.1 to 1.3:         -   1.1 Moving the position indicator in a fixed position.         -   1.2 Detecting the sensor signals a₁ to a_(r) and storing as             vector {right arrow over (a)}_(F).         -   1.3 Executing the following steps for a predetermined number             v=1 to b≧1 of different external fields, wherein the number             b fulfills the condition b≦r−2:             -   1.3.1 Applying the external field.                 -   The external field is for example a homogeneous                     magnetic field, which points in an arbitrary space                     direction. The external field can also be an                     inhomogeneous external field, to which the position                     measuring system is exposed. In order to eliminate                     the influence of external fields, additional sensors                     4 are required. The minimum number of sensors 4,                     which is necessary to eliminate a number of b                     different external fields, is according to the above                     equation r=b+2.             -   1.3.2 Detecting the sensor signals a₁ to a_(r) and                 forming an external field vector {right arrow over                 (ƒ)}_(v)={right arrow over (a)}−{right arrow over                 (a)}_(F)                 -   or given the case of an offset compensated external                     field vector, respectively, {right arrow over                     (ƒ)}_(v)={right arrow over (a)}−{right arrow over                     (a)}_(F)−{right arrow over (a)}₀,             -   In the external field vector {right arrow over (ƒ)}_(v),                 the proportion of the sensor signals generated by the                 position indicator is thus eliminated.

Afterwards, the steps according to variant 1 or 2 follow:

In variant 1, the dimension of the vector space of the measurement vectors is reduced. Variant 1 comprises the steps:

-   -   2. Forming a complete, orthonormal basis of the measurement         vector space with the steps         -   2.1 Forming of c orthonormal vectors {right arrow over             (g)}₁, . . . , {right arrow over (g)}_(c) of the vector             space spanned by the external field vectors {right arrow             over (ƒ)}₁, . . . , {right arrow over (ƒ)}_(b), wherein c≦b.             Thus, for these basis vectors is {right arrow over             (g)}_(i)·{right arrow over (g)}_(j)=0 i≠j and {right arrow             over (g)}_(i)·{right arrow over (g)}_(i)=1 for i=1, . . . ,             c.             -   The orthonormal basis of the vectors {right arrow over                 (g)}₁, . . . , {right arrow over (g)}_(c) can for                 example be calculated with the                 Gram-Schmidt-orthogonalisation method, which is a                 standard method of linear algebra and which                 orthonormalizes the external field vectors {right arrow                 over (ƒ)}₁, . . . , {right arrow over (ƒ)}_(b)                 step-by-step.         -   2.2 Complementing the c vectors {right arrow over (g)}₁, . .             . , {right arrow over (g)}_(c) with n=r−c further vectors             {right arrow over (g)}_(c+1), . . . , {right arrow over             (g)}_(c+n) to a complete, orthonormal basis of the             r-dimensional vector space of the signal vectors. Thus, for             these basis vectors is {right arrow over (g)}_(i)·{right             arrow over (g)}_(j)=0 for i≠j and {right arrow over             (g)}_(i)·{right arrow over (g)}_(i)=1 for i=1, . . . , r.     -   3. Forming the projection matrix P according to

$P = {\begin{pmatrix} {\overset{\rightarrow}{g}}_{c + 1} \\ \vdots \\ {\overset{\rightarrow}{g}}_{c + n} \end{pmatrix} = \begin{pmatrix} g_{{c + 1},1} & \ldots & g_{{c + 1},r} \\ \vdots & \; & \vdots \\ g_{{c + n},1} & \ldots & g_{{c + n},r} \end{pmatrix}}$

-   -   -   The n basis vectors {right arrow over (g)}_(c+1), . . . ,             {right arrow over (g)}_(c+n) form the rows of the n×r             matrix P. The application of P to a signal vector {right             arrow over (a)} delivers a n-dimensional vector {right arrow             over (q)}=P{right arrow over (a)} with respect to the new             basis {right arrow over (g)}_(c+k), k=1, . . . , n, which             contains no portions of external field vectors. In this             case, when calculating the measurement vectors the dimension             r of the vector space of the signal vectors is reduced by             the value c, so that the space of the measurement vectors             has the dimension n=r−c.

In variant 2, the dimension of the vector space of the measurement vectors is not reduced compared to the dimension of the vector space of the signal vectors. The variant 2 comprises the steps:

-   -   2. Forming an orthonormal basis of c vectors {right arrow over         (g)}₁, . . . , {right arrow over (g)}_(c) of the vector space         spanned by the external field vectors {right arrow over (ƒ)}₁, .         . . , {right arrow over (ƒ)}_(b), wherein c≦b. Therefore, for         these basis vectors it is {right arrow over (g)}_(i)·{right         arrow over (g)}_(j)=0 for i≠j and {right arrow over         (g)}_(i)·{right arrow over (g)}_(i)=1 for i=1, . . . , c.         -   The orthonormal basis of the vectors {right arrow over             (g)}₁, . . . , {right arrow over (g)}_(c) can for example be             calculated with the Gram-Schmidt-orthogonalisation method,             which is a standard method of linear algebra and which             orthonormalizes the external field vectors {right arrow over             (ƒ)}₁, . . . , {right arrow over (ƒ)}_(b) step-by-step.     -   3. Forming the projection matrix P according to

${P = { - {\sum\limits_{i = 1}^{c}{{\overset{\rightarrow}{g}}_{i} \otimes {\overset{\rightarrow}{g}}_{i}}}}},$

-   -   -   wherein {right arrow over (g)}_(i)             {right arrow over (g)}_(i) are the tensor products of the             vectors {right arrow over (g)}_(i), i.e. are the r×r             matrices

$\mspace{20mu} {{{\overset{\rightarrow}{g}}_{i} \otimes {\overset{\rightarrow}{g}}_{i}} = {{{\begin{pmatrix} {g_{i\; 1}g_{i\; 1}} & \ldots & {g_{i\; 1}g_{ir}} \\ \vdots & \; & \vdots \\ {g_{ir}g_{i\; 1}} & \ldots & {g_{ir}g_{ir}} \end{pmatrix}.\mspace{20mu} {The}}\mspace{14mu} {equation}\mspace{14mu} P} = { - {\sum\limits_{i = 1}^{c}\; {{{\overset{\rightarrow}{g}}_{i} \otimes {\overset{\rightarrow}{g}}_{i}}\mspace{14mu} {reads}\mspace{14mu} {written}\mspace{14mu} {in}\mspace{14mu} {full}\text{:}}}}}}$ $\begin{pmatrix} P_{11} & \ldots & P_{1r} \\ \vdots & \; & \vdots \\ P_{r\; 1} & \ldots & P_{rr} \end{pmatrix} = {\begin{pmatrix} 1 & \; & \; \\ \; & \ddots & \; \\ \; & \; & 1 \end{pmatrix} - \begin{pmatrix} {g_{11}g_{11}} & \ldots & {g_{11}g_{1r}} \\ \vdots & \; & \vdots \\ {g_{1r}g_{11}} & \ldots & {g_{1r}g_{1r}} \end{pmatrix} - \ldots - \begin{pmatrix} {g_{c\; 1}g_{c\; 1}} & \ldots & {g_{c\; 1}g_{cr}} \\ \vdots & \; & \vdots \\ {g_{cr}g_{c\; 1}} & \ldots & {g_{cr}g_{cr}} \end{pmatrix}}$

-   -   -   The application of the protection matrix P on a signal             vector {right arrow over (a)} delivers an r-dimensional             vektor {right arrow over (q)}=P{right arrow over (a)} with             respect to the standard base of the measurement vector             space, wherein the vector {right arrow over (q)} is             orthogonal to the subspace of the external field vectors.             The dimension of the vector space of the measurement vector             {right arrow over (q)} is equal to the dimension of the             vector space of the signal vectors {right arrow over (a)}.             Therefore, the projection matrix P is a r×r matrix.

Determination of the Matrix M

The determination of the matrix M occurs in a calibration phase. For this, two methods are explained in the following. The first method is based on the condition that it is ensured “from outside” that calibration positions are selected, which are distributed over the whole measurement range and enable a precise result. The second method is a bootstrap-process, which is started with a factory setting and then iteratively produces better calibration positions.

Method 1

The determination of the matrix M occurs with the following steps:

-   -   1. Detecting the sensor signals a_(k1) to a_(kr) of the r         sensors 4 for a predetermined number of k=1 to m positions of         the position indicator, wherein m≧3.         -   This step delivers in signal vectors {right arrow over             (a)}₁, . . . , {right arrow over (a)}_(m) with {right arrow             over (a)}_(k)=(a_(k1), . . . , a_(kr)), k=1 to m.     -   2. Forming nm measurement vectors {right arrow over (q)}_(k)         with k=1 to m according to the predetermined function g to         {right arrow over (g)}_(k)={right arrow over (g)}({right arrow         over (a)}_(k)),     -   3. Forming a matrix Q according to

$Q = {\begin{pmatrix} {\overset{\rightarrow}{q}}_{1} \\ \vdots \\ {\overset{\rightarrow}{q}}_{m} \end{pmatrix} = {\begin{pmatrix} q_{11} & \ldots & q_{1n} \\ \vdots & \; & \vdots \\ q_{m\; 1} & \ldots & q_{mn} \end{pmatrix}.}}$

-   -   -   The in measurement vectors {right arrow over (q)}₁, . . . ,             {right arrow over (q)}_(m) form the rows of the matrix Q.             The matrix Q is a m×n matrix, wherein the number n of its             columns depends on the predetermined function g. It is n≦r.             As explained further above, a reduction to n<r can occur,             for example, during the elimination of external field             influences.         -   The determination of the matrix M from the matrix Q occurs             in such a way that the transformation {right arrow over             (p)}=M·{right arrow over (q)} is a linear mapping, which             maps each of the m measurement vectors {right arrow over             (q)}_(k) on a 2-component vector {right arrow over (p)}_(k),             wherein the tips of the vectors {right arrow over (p)}_(k)             substantially lie on a circle or on a circular arc.             “Substantially” means in ideal circumstances really, in             practice, however, usually only approximately, i.e. the tips             {right arrow over (p)}_(k) of the vectors are then             distributed on both sides of the circle line. The circle can             therefore be a dented circle, with bumps inside and outside.             Because the transformation {right arrow over (p)}=M·{right             arrow over (q)} is a linear mapping, the matrix M can be             determined with methods of linear algebra, wherein these             methods include a decomposition of the matrix Q in two or             three matrices. The determination of the matrix M requires             no statistical methods as for example regression analysis or             fits. The determination of the matrix M from the matrix Q             comprises the following steps 4 to 9:

    -   4. Determining three matrices U, S and V^(T), which fulfill the         equations Q=U·S·V^(T), wherein U and V are orthogonal matrices,         i.e. fulfill the equations

U ^(T) ·U=

V ^(T) ·V=

, and

-   -   -   S is a positive semidefinite m×n diagonal matrix. U^(T) and             V^(T) are the transposed matrices of the matrices U and V,             respectively.             designates the unit matrix. The equation Q=U·S·V^(T) reads             written in full:

$\begin{pmatrix} q_{11} & \ldots & q_{1n} \\ \vdots & \; & \vdots \\ q_{m\; 1} & \ldots & q_{mn} \end{pmatrix} = {\begin{pmatrix} U_{11} & \ldots & U_{1m} \\ \vdots & \; & \vdots \\ U_{m\; 1} & \ldots & U_{mm} \end{pmatrix} \cdot \begin{pmatrix} \begin{matrix} s_{1} & \; & \; \\ \; & \ddots & \; \\ \; & \; & s_{l} \end{matrix} & \begin{matrix} 0 & \ldots & 0 \\ \vdots & \; & \vdots \\ 0 & \ldots & 0 \end{matrix} \\ \begin{matrix} 0 & \ldots & 0 \\ \vdots & \; & \vdots \\ 0 & \ldots & 0 \end{matrix} & \begin{matrix} 0 & \ldots & 0 \\ \vdots & \; & \vdots \\ 0 & \ldots & 0 \end{matrix} \end{pmatrix} \cdot \begin{pmatrix} V_{11}^{T} & \ldots & V_{1n}^{T} \\ \vdots & \; & \vdots \\ V_{n\; 1}^{T} & \ldots & V_{nn}^{T} \end{pmatrix}}$

-   -   -   The matrix U is a m×m matrix, the matrix V is a n×n matrix.             The values s₁ to s_(l), l≦min(m,n) of the diagonal values of             the matrix S are called singular values of the matrix i.e.             s_(j)≧0 for j=1, . . . , l≦min(m,n). (V_(jk) ^(T) designates             the jk-element of the transposed matrix V^(T))         -   The determination of the matrix M does not require all, but             only selected and non-zero elements of the diagonal matrix             S, as well as elements of the matrices V and U linked to             these. The matrices U, S and V are therefore advantageously             calculated with the aid of the singular value decomposition             (SVD), a standard operation of linear algebra, because the             singular value decomposition sorts the elements of the             diagonal matrix S according to their size, so that s₁≧s₂≧ .             . . ≧s_(l). The singular value decomposition delivers the             matrix V, the matrix V^(T) is obtained by transposition.         -   The determination of the matrix M is performed as a function             of an integer parameter h, which must fulfill the condition             2≦h≦l. The accuracy of the matrix M is the higher, the             higher the value of the parameter h is.

    -   5. Forming sub-matrices Ũ, {tilde over (S)} and {tilde over         (V)}^(T), wherein         -   Ũ is that m×h sub-matrix of the m×m matrix U, which is             formed by discarding the columns h+1 to m,         -   {tilde over (S)} is that h×h sub-matrix of the m×n matrix S,             which is formed by discarding the rows h+1 to m and             discarding the columns h+1 to n, and         -   {tilde over (V)}^(T) is that h×n sub-matrix of the n×n             matrix V^(T), which is formed by discarding the rows h+1 to             n.

    -   6. Forming a 2×m auxiliary matrix W:         -   If the measurement vectors are present in an arbitrary             order, then the auxiliary matrix W is formed according to             variant 1. If the measurement vectors {right arrow over             (q)}_(k) with k=1, . . . , m are sorted in ascending or             descending order according to their associated position             φ_(k), then the auxiliary matrix W is formed according to             variant 2.         -   Variant 1 comprises the sub steps:         -   6.1 For i=1 to m calculating α_(i)=a tan             2(U₁₂U_(i1)−U₁₁U_(i2), U₁₂U_(i2)+U₁₁U_(i1)).         -   6.2 Calculating the permutation k(i), which sorts the             α_(i)-values according to their size, so that             α_(k(i))>α_(k(i−1)) for i=2, . . . , m.         -   6.3 Forming the auxiliary matrix W according to

$W_{i,{k{(j)}}} = \left\{ {\begin{matrix} {{\sin \left( \frac{2\; {\pi \left( {j - 1} \right)}}{m} \right)},{i = 1}} \\ {{\cos \left( \frac{2\; {\pi \left( {j - 1} \right)}}{m} \right)},{i = 2}} \end{matrix},{{{for}\mspace{14mu} j} = {1\mspace{14mu} {to}\mspace{14mu} m}},} \right.$

-   -   -    wherein W_(i,k(j)) designates the i,k(j)-element of the             matrix W, i.e. in the column k(j) given by the previously             determined permutation of the indices is the entry for the             index j. In other words, the columns of the matrix W are             sorted depending on the order of the measurement vectors.         -   Variant 2 comprises the single step depending on the             auxiliary parameter ƒ=1 for ascending and ƒ=−1 for             descending order according to the associated position φ_(k)         -   6.1 Forming the auxiliary matrix W according to

$W_{ij} = \left\{ {\begin{matrix} {{\sin \left( {f \cdot \frac{2\; {\pi \left( {j - 1} \right)}}{m}} \right)},{i = 1}} \\ {{\cos \left( {f \cdot \frac{2\; {\pi \left( {j - 1} \right)}}{m}} \right)},{i = 2}} \end{matrix},{{{for}\mspace{14mu} j} = 1}} \right.$

-   -   -    to m

    -   7. Calculating a 2×h auxiliary matrix T according to T=WŨ         -   This equation reads written in full:

$\begin{pmatrix} T_{11} & \ldots & T_{1h} \\ T_{21} & \ldots & T_{2\; h} \end{pmatrix} = {\begin{pmatrix} W_{11} & \ldots & W_{1m} \\ W_{21} & \ldots & W_{2\; m} \end{pmatrix} \cdot \begin{pmatrix} U_{11} & \ldots & U_{1h} \\ \vdots & \; & \vdots \\ U_{m\; 1} & \ldots & U_{mh} \end{pmatrix}}$

-   -   8. Forming a 2×n auxiliary matrix H         -   The auxiliary matrix H is formed from the auxiliary matrix T             and the sub-matrices {tilde over (S)} and {tilde over             (V)}^(T) according to H=T{tilde over (S)}⁻¹{tilde over             (V)}^(T). This equation reads written in full:

$\begin{pmatrix} H_{11} & \ldots & H_{1n} \\ H_{21} & \ldots & H_{2\; n} \end{pmatrix} = {\begin{pmatrix} T_{11} & \ldots & T_{1h} \\ T_{21} & \ldots & T_{2\; h} \end{pmatrix} \cdot \begin{pmatrix} \frac{1}{s_{1}} & \; & \; \\ \; & \ddots & \; \\ \; & \; & \frac{1}{s_{h}} \end{pmatrix} \cdot \begin{pmatrix} V_{11}^{T} & \ldots & V_{1n}^{T} \\ \vdots & \; & \vdots \\ V_{h\; 1}^{T} & \ldots & V_{hn}^{T} \end{pmatrix}}$

-   -   9. Forming the matrix M         -   With the calculation of the auxiliary matrix H all steps for             forming the matrix M are completed, with the exception of a             single step, which, however, is not required in all cases.             This step concerns the determination or check of the             movement direction. In an angle sensor, the movement             direction corresponds to the sense of rotation. Note the             following:         -   If in normal operation only the absolute value |φ₁−φ₂| of             the difference of two measured positions φ₁ and φ₂ is             important, then it is M=H.         -   If in normal operation also the movement direction, i.e. the             sign of the difference φ₁−φ₂ is important, then it is M=NH,             where

$N = \begin{pmatrix} d & 0 \\ 0 & 1 \end{pmatrix}$

-   -   -    is an auxiliary matrix, whose parameter d has the value d=1             or d=−1.         -   If the measurement vectors are sorted according to ascending             or descending positions, it is d=1, or N=             , respectively (             denotes the unit matrix) and therefore M=H. The correct             movement direction results here therefrom, that the             distinction between ascending and descending positions             enters into the calculation of the auxiliary matrix W by the             parameter ƒ. Therefore, the auxiliary matrix W contains the             information on the movement direction.         -   If the measurement vectors are not sorted according to             ascending or descending positions, but are provided in an             arbitrary order, then the movement direction can be             determined by a measurement with the substeps 9.1 to 9.4:         -   9.1 Bringing the position indicator in a first predetermined             set-position ψ₁, detecting the sensor signals and forming             the corresponding measurement vector, and calculating the             corresponding actual position φ₁, and bringing the position             indicator in a second predetermined set-position ψ₂,             detecting the sensor signals and forming the corresponding             measurement vector, and calculating the corresponding actual             position φ₂, wherein each time the calculation occurs with             the steps C and D of the normal operation, with use of the             matrix M=H.         -   9.2 Calculating the difference ψ₁−ψ₂,         -   9.3 Calculating the difference φ₁−φ₂, and         -   9.4 Setting the parameter d=+1 if the sign of the difference             φ₁−φ₂ equals the sign of the difference ψ₁−ψ₂, and d=−1             otherwise.         -   The formation of the matrix M can formally always be written             as a M=NH, because the auxiliary matrix N can be the unit             matrix             . The movement direction or the parameter d characterizing             it can also always be determined by a measurement with the             substeps 9.1 to 9.4. The auxiliary matrix N leaves the             movement direction as it is or reverses it.

Special Cases

The described method can also be used with minor modifications in the following special cases:

-   -   1. Angle measuring system for a predetermined angular range         -   The positions are in a limited interval φ∈[φ_(min),φ_(max)]             instead of φ∈[0,2π]. (For a linear position meter, this             corresponds to the use of only a fraction of the full period             distance D).         -   When the measurement vectors are disordered, the elements of             the auxiliary matrix W are calculated according to

$W_{i,{k{(j)}}} = \left\{ {\begin{matrix} {{\sin \left( \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m} \right)},{i = 1}} \\ {{\cos \left( \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m} \right)},{i = 2}} \end{matrix},{{{for}\mspace{14mu} j} = {1\mspace{14mu} {to}\mspace{14mu} m}},} \right.$

-   -   -   wherein φ_(min) and φ_(max) denote the lower or upper limit             of the angular range used, and k(j) represents the             permutation, which orders the measurement vectors {right             arrow over (q)}_(j) in ascending or descending order             according to their position,         -   and when the measurement vectors are ordered depending on             the auxiliary quantity ƒ=1 for ascending and ƒ=−1 for             descending order according to

$W_{ij} = \left\{ {\begin{matrix} {{\sin \left( {f \cdot \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m}} \right)},{i = 1}} \\ {{\cos \left( {f \cdot \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m}} \right)},{i = 2}} \end{matrix},{{{for}\mspace{14mu} j} = {1\mspace{14mu} {to}\mspace{14mu} {m.}}}} \right.$

-   -   2. Angle measurement system with a position indicator which has         a multipole character         -   The position sensor can have a multipole character, in             particular if the position indicator is a magnet, but also             when the position indicator has coils. The multipole             character of the magnet or of the coils is characterized by             a parameter, e, which indicates the number of repetitions of             the field distribution during a complete rotation of the             position indicator. For a magnet with a dipolar field it is             e=1, while for a quadrupole field it is e=2.         -   The elements of the auxiliary matrix W are calculated in the             case of ordered measurement vectors as a function of the             auxiliary quantity f=1 for ascending and f=−1 for descending             order according to

$W_{ij} = \left\{ {\begin{matrix} {{\sin \left( {{fe} \cdot \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m}} \right)},{i = 1}} \\ {{\cos \left( {{fe} \cdot \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m}} \right)},{i = 2}} \end{matrix},{{{for}\mspace{14mu} j} = {1\mspace{14mu} {to}\mspace{14mu} m}},} \right.$

-   -   -   and for disordered measurement vectors according to

$W_{i,{k{(j)}}} = \left\{ {\begin{matrix} {{\sin \left( {e \cdot \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m}} \right)},{i = 1}} \\ {{\cos \left( {e \cdot \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m}} \right)},{i = 2}} \end{matrix},{{{for}\mspace{14mu} j} = {1\mspace{14mu} {to}\mspace{14mu} m}},} \right.$

-   -   -   The calculation of the position φ from the 2-component             vector {right arrow over (p)} occurs in this case according             to

$\varphi = {\frac{f\left( \overset{\rightarrow}{p} \right)}{e} = {\frac{f\left( {p_{1},p_{2}} \right)}{e}.}}$

Determination of the Matrix M in the Case n=2

The determination of the matrix M for n=2 can be done in three different ways, namely according to the general method with the previously described steps 1-9, or in a simplified manner according to one of the following two methods:

Method 1 comprises:

-   -   the steps 1 to 5 of the general method,     -   forming an auxiliary matrix H as H={tilde over (S)}⁻¹{tilde over         (V)}^(T), and     -   forming the matrix M according to variant 1 or variant 2:         -   Variant 1         -   If in normal operation only the absolute value |φ₁−φ₂| of             the difference of two measured positions φ₁ and φ₂ is             important, then it is M=H         -   Variant 2         -   If in normal operation also the movement direction, i.e. the             sign of the difference φ₁−φ₂ is important, then it is M=NH.         -   If the measurement vectors are sorted according to ascending             or descending positions, the parameter d can be calculated             with the auxiliary quantity f=1 for ascending or f=−1 for             descending order and from selected elements of the matrix U             to

$d = {f \cdot {\frac{{U_{12}U_{21}} - {U_{11}U_{22}}}{{{U_{12}U_{21}} - {U_{11}U_{22}}}}.}}$

-   -   -   If the measurement vectors are not sorted according to             ascending or descending positions but are in any order, then             the parameter d can be determined by the above-described             measurement with the substeps 9.1 to 9.4.

With method 1 there is no need to determine the auxiliary matrices W and T.

Method 2 comprises the steps 1 to 3 of the general method. This delivers the m×2 Matrix Q. The steps 4 to 8 of the general method are replaced by the following two steps 4 and 5:

-   -   4. Determining two matrices U and B, which fulfill the equation         Q=U·B, wherein U is a m×2 matrix with U^(T)·U=         , and B is a symmetrical, positive definite 2×2 matrix. The         equation Q=U·B reads written in full:

$\begin{pmatrix} q_{11} & \ldots & q_{12} \\ \vdots & \; & \vdots \\ q_{m\; 1} & \ldots & q_{m\; 2} \end{pmatrix} = {\begin{pmatrix} U_{11} & \ldots & U_{12} \\ \vdots & \; & \vdots \\ U_{m\; 1} & \ldots & U_{m\; 2} \end{pmatrix} \cdot \begin{pmatrix} B_{11} & B_{12} \\ B_{21} & B_{22} \end{pmatrix}}$

-   -   -   The matrix B may be calculated by the substeps:         -   4.1 Forming a symmetrical, positive definite 2×2 auxiliary             matrix D to D=Q^(T)·Q,         -   4.2 Calculating the matrix B as square root according to             B=D^(1/2). The calculation of the square root of a matrix is             a standard procedure of linear algebra.         -   The matrix U results to U=QB⁻¹, wherein the matrix U is             necessary, however, at most for the calculation of the             movement direction.

    -   5. Forming an auxiliary matrix H according to H=B⁻¹.

This is followed by step 9 of the general method, with the modification described below, when the measurement vectors are arranged in ascending or descending order according to their positions. Because this method does not calculate the auxiliary matrix W, the correct movement direction does not enter automatically. However, it can be calculated with the auxiliary quantity ƒ=1 for ascending or ƒ=−1 for descending order and from selected elements of the matrix U to:

$d = {f \cdot {\frac{{U_{12}U_{21}} - {U_{11}U_{22}}}{{{U_{12}U_{21}} - {U_{11}U_{22}}}}.}}$

Absolute Positions

If, in normal operation, absolute positions are also to be determined, a determination of the zero position is necessary, for example by the following steps:

-   -   positioning the position indicator in the zero position,     -   determining the position φ₀ with the steps A to D of the normal         operation, and     -   changing the function ƒ({right arrow over (p)}) to ƒ({right         arrow over (p)})=ƒ({right arrow over (p)})−φ₀.

FIGS. 6 to 8 illustrate the described method steps. On the left side of the FIGS. 6 to 8 and also in FIGS. 9 and 10 explained further below, the position measuring system is represented symbolically. In detail:

FIG. 6 illustrates the normal operation of the position measuring system. The position indicator is in a position φ and generates the sensor signals a₁, . . . , a_(r) in the sensors 4. These signals are offset compensated by the offset vector {right arrow over (a)}₀ and external field compensated by the projection matrix P (each optionally), which delivers the measurement vector {right arrow over (q)}. Subsequently, the 2-component vector {right arrow over (p)} is calculated by multiplying the matrix M with the measurement vector {right arrow over (q)}. Finally from the components of {right arrow over (p)} the position φ is calculated according to φ=a tan 2(p₂,p₁), or according to another function.

FIG. 7 illustrates the determination of the projection matrix P in the calibration phase. For doing so, b external fields are applied to the position measuring system, which deliver the signal vectors {right arrow over (a)}_(j)=(a_(j1), . . . , a_(jr)) j=1, . . . , b for the external fields. From the signal vectors {right arrow over (a)}_(j), optionally, offset compensated external field vectors {right arrow over (ƒ)}_(j) are formed. From this, the orthonormal basis {right arrow over (g)}₁, . . . , {right arrow over (g)}_(c) of the subspace spanned by the external field vectors is formed, for example by the Gram-Schmidt method GS.

The projection matrix P is calculated from the basis vectors {right arrow over (g)}_(j), in that the matrices {right arrow over (g)}_(j)

{right arrow over (g)}_(j) obtained as tensor products from the basis vectors are subtracted from the unit matrix.

FIG. 8 illustrates the determination of the matrix M in the calibration phase. For this, the signal vectors {right arrow over (a)}_(k)=(a_(k1), . . . , a_(kr)), k=1, . . . , m are compensated with the offset vector {right arrow over (a)}₀ and the projection matrix P (each optionally). The resulting measurement vectors {right arrow over (q)}_(k) form the rows of the matrix Q. The matrix Q is then decomposed into the three matrices Ũ, {tilde over (S)} and {tilde over (V)}^(T). Depending on the predetermined parameter h, the sub-matrices Ũ, {tilde over (V)}^(T) and {tilde over (S)} are determined from U, V and S. The matrix M is then calculated with the auxiliary matrices T and N. For this purpose, the matrix T=WŨ is formed from the auxiliary matrices W and Ũ.

Method 2: Bootstrap Method

The determination of the matrix M can also take place by means of a bootstrap method. The bootstrap process includes an initialization on the manufacturer's side, as well as a measurement data acquisition during commissioning and an iterative process for increasing the accuracy of the matrix M.

The bootstrap method comprises the following steps:

-   -   1. Setting an index i=0 and configuring a predetermined initial         matrix M⁽⁰⁾.         -   The initial matrix M⁽⁰⁾ is a manufacturer-side             initialization of the matrix M=M⁽⁰⁾, so that—possibly             inaccurate—position calculations can be carried out with the             position measuring system.     -   2. Detecting and storing a plurality of measurement vectors         {right arrow over (q)} according to the steps A and B of the         normal operation.     -   3. Iterative determination of matrices M^((i)) (with i>0) by         repeatedly performing the following steps:         -   3.1. setting the matrix M to M=M^((i))         -   3.2. increasing the index i by 1         -   3.3. calculating the positions φ of the stored measurement             vectors {right arrow over (q)} according to the steps C and             D of the normal operation,         -   3.4. selecting a number of k=1 to nm of the stored             measurement vectors {right arrow over (q)}_(k), the             corresponding positions φ_(k) of which fulfill a             predetermined criterion,         -   3.5. forming a matrix Q according to

$Q = \begin{pmatrix} {\overset{->}{q}}_{1} \\ \vdots \\ {\overset{->}{q}}_{m} \end{pmatrix}$

-   -   -   3.6. determining a new matrix M^((i)) from the matrix Q in             such a way, that the transformation {right arrow over             (q)}=M·{right arrow over (q)} is a linear mapping, which             maps each of the m selected measurement vectors {right arrow             over (q)}_(k) on a 2-component vector {right arrow over             (p)}_(k), until the matrix M^((i)) within a predetermined             tolerance no longer differs from the matrix M^((i−1)) of the             previous iteration, wherein the tips of the vectors {right             arrow over (p)}_(k) with increasing index i substantially             come to lie on a circle 6 or on a circular arc. This             terminates the bootstrap process (either with M=M^((i)) or             M=M^((i-1))).         -   Subsequently, a setting or calibration of the zero position             follows, as described above.

The individual iterations adaptively lead to improved matrices M^((i)), because in each step more precise positions φ of the intermediately stored measurement vectors {right arrow over (q)} can be determined, so that the matrix Q is formed according to the predetermined criterion from more suitable measurement vectors {right arrow over (q)} and thus leads to a more precise matrix M. The desired sense of rotation is already included in the start matrix M and does not change during the iterations.

Recalibration

The method also allows recalibration during normal operation, which is carried out in two steps (and is similar to the bootstrap process):

-   -   1. Capturing measurement results for the recalibration with the         steps A to D of the normal operation with aid of the substeps:         -   1.1 Capturing and buffering measurement vectors {right arrow             over (q)} according to the steps A and B for a plurality of             measurements, and         -   1.2 Calculating the corresponding positions φ according to             the steps C and D and buffering these positions φ for all             buffered measurement vectors {right arrow over (q)}.     -   2. The actual recalibration takes place with the following steps         2.1 to 2.4         -   2.1 Selecting m buffered measurement vectors {right arrow             over (q)}, the corresponding positions φ of which fulfill a             predetermined criterion,         -   2.2 Rebuilding the matrix Q from the selected measurement             vectors {right arrow over (q)}.         -   2.3 Rebuilding the matrix M, from the matrix Q, and         -   2.4 given the case, redetermining the zero position.

The predetermined criterion comprises preferably one or more of the following criteria:

-   -   the positions φ are distributed over the entire range from 0 to         2π,     -   the positions φ are equally distributed within a predetermined         tolerance criterion,     -   the measurement vectors {right arrow over (q)} are sorted in         ascending or in descending order of the positions φ.

A large number of measurement vectors improves the calibration when the positions φ are not equidistributed.

Uniform distribution within a given tolerance criterion means that the k=1 to m positions φ_(k) of the selected measurement vectors {right arrow over (q)}_(k) correspond up to a defined tolerance Δφ to the values

$\frac{{\left( {k - 1} \right) \cdot 2}\pi}{m},$

i.e. it is

${{\varphi_{k} - \frac{{\left( {k - 1} \right) \cdot 2}\pi}{m}}} < {{\Delta\varphi}.}$

Plausibility Check

The position of the sensor chip 3 may change over time relative to the position of the rotation axis 1 or the axis x. In order to recognize this, the method can be extended by a plausibility check, as long as the dimension n of the vector space of the measurement vectors is n>2. The method requires a further n×n matrix C, which is determined in the calibration phase and used in normal operation for the plausibility check.

The calibration phase then additionally comprises the following step:

-   -   Forming the matrix C from the sub-matrix {right arrow over (V)}         according to

C=

−{tilde over (V)}·{tilde over (V)} ^(T)

-   -   This equation reads written in full:

$\begin{pmatrix} C_{11} & \ldots & C_{1n} \\ \vdots & \; & \vdots \\ C_{n\; 1} & \ldots & C_{nn} \end{pmatrix} = {\begin{pmatrix} 1 & \; & \; \\ \; & \ddots & \; \\ \; & \; & 1 \end{pmatrix} - {\begin{pmatrix} V_{11} & \ldots & V_{1h} \\ \vdots & \; & \vdots \\ V_{n\; 1} & \ldots & V_{nh} \end{pmatrix} \cdot \begin{pmatrix} V_{11}^{\top} & \ldots & V_{1n}^{\top} \\ \vdots & \; & \vdots \\ V_{h\; 1}^{\top} & \ldots & V_{hn}^{\top} \end{pmatrix}}}$

The normal operation additionally comprises the following step E:

-   -   E) calculating a check value with the substeps E.1 to E.3:         -   E.1 calculating a vector {right arrow over (d)} to

{right arrow over (d)}=C{right arrow over (q)}

-   -   -   This equation reads written in full:

$\begin{pmatrix} d_{1} \\ \vdots \\ d_{n} \end{pmatrix} = {\begin{pmatrix} C_{11} & \ldots & C_{1n} \\ \vdots & \; & \vdots \\ C_{n\; 1} & \ldots & C_{nn} \end{pmatrix} \cdot \begin{pmatrix} q_{1} \\ \vdots \\ q_{n} \end{pmatrix}}$

-   -   -   E.2 calculating a check value v from the vector {right arrow             over (d)}             -   The check value can for example be the quantity

$v = {\sqrt{\frac{\sum\limits_{l = 1}^{N_{s}}d_{l}^{2}}{\sum\limits_{l = 1}^{N_{s}}q_{l}^{2}}}.}$

-   -   -   -    The check value can also be another quantity, in                 particular one which has the properties of a vector                 norm. A vector norm typically has the following                 properties:             -   (i) v({right arrow over (d)})=0 exactly then, when                 {right arrow over (d)}={right arrow over (0)} (zero                 vector),             -   (ii) for each R∈                 (set of real numbers) is v(R·{right arrow over                 (d)})=|R|·v({right arrow over (d)}), and             -   (iii) for all vectors {right arrow over (d)}₁,{right                 arrow over (d)}₂ is v({right arrow over (d)}₁+{right                 arrow over (d)}₂)≦v({right arrow over (d)}₁)+v({right                 arrow over (d)}₂).                 -   The check value can also be a quantity, which only                     fulfills the properties (i) and (ii) of such a                     vector norm.

        -   E.3 triggering a predetermined action, if the check value v             exceeds a predetermined value.

FIG. 9 illustrates the determination of the matrix C in the calibration phase. This calculation proceeds until the determination of the three matrices Ũ, {tilde over (S)} and {tilde over (V)}^(T) analogously to the determination of the matrix M. The matrix C is then calculated, in that the product of the sub-matrix {tilde over (V)} of the matrix V with its transposed {tilde over (V)}^(T) is subtracted from the unit matrix.

FIG. 10 illustrates the plausibility check. This check is carried out in normal operation. This calculation proceeds until the determination of {right arrow over (q)} analogously to the calculation of the position (FIG. 6). For the plausibility check, the n-component vector {right arrow over (d)} is calculated from the measurement vector {right arrow over (q)} by means of the matrix C and from this finally the check value v.

Position Determination and External Field Suppression in the Special Case n=2 and r>n

The determination of the position φ of a rotating or linearly moving position indicator may include an external field compensation, with which the projection matrix P maps the signal vectors of the dimension r on measurement vectors of the dimension n=2, and a calculation of the position a from the components q₁ and q₂ of a measurement vector {right arrow over (q)} according to a classical method. The determination of the projection matrix P takes place in a calibration phase according to the above-described variant 1 with the values b=r−2 and c=b.

In normal operation the determination of the position φ takes place as follows:

-   -   A) Detecting the sensor signals a₁ to a_(r) of the sensors (4)         and forming a signal vector {right arrow over (a)}=((a₁, . . . ,         a_(r)).     -   B) Forming a measurement vector {right arrow over (q)} according         to {right arrow over (q)}=g({right arrow over (a)}), wherein         -   g({right arrow over (a)})=P·{right arrow over (a)}, if             (only) an external field correction is performed, or         -   g({right arrow over (a)})=P·({right arrow over (a)}−{right             arrow over (a)}₀), if an offset and an external field             correction is performed.         -   The vector {right arrow over (q)} is a 2-component vector             with the components q₁ and q₂. The calculation of the vector             {right arrow over (p)} according to {right arrow over             (p)}=M·{right arrow over (q)} simplifies in this case to             {right arrow over (p)}={right arrow over (q)}, because in             this special case the matrix M is the unit matrix.     -   If the movement of the position indicator provides a sequence of         measurement vectors whose tips lie on a circle or circular arc,         then the step follows:     -   C) Determining the position φ from the component p₁ and p₂ of         the vector {right arrow over (p)} by means of a predetermined         function ƒ({right arrow over (p)}) to φ=ƒ({right arrow over         (p)})=ƒ(p₁,p₂), wherein the predetermined function ƒ(q=ƒ(q₁,q₂)         is based on the equations q₁=cos(φ) and q₂=sin(φ).     -   If the movement of the position indicator provides a sequence of         measurement vectors whose tips lie on an ellipse, then the         calibration phase has to be supplemented such that the         parameters of the ellipse are determined after the determination         of the projection matrix P, so that then the calculation of the         position φ can take place on the basis of the values q₁ and q₂         and the parameters of the ellipse. The measurement value pairs         q₁ and q₂ of the measurement vector {right arrow over (q)} can         be processed in the calibration phase as well as in normal         operation, for example, according to the method of the         international publication WO 2005124286 A1, wherein the values         q₁ and q₂ correspond to the measurement values which are         designated there as measurement value pairs (x_(i),y_(i)) with         i=1, . . . , N, or according to any other known method for         approximately sine or cosine shaped measurement signals. I.e.         also here a function ƒ({right arrow over (q)}) can be determined         which delivers the position φ to φ=ƒ({right arrow over (q)}).

This method enables the separation of external field detection and position determination in a manner in which it is not necessary to position the sensors very precisely and to adjust their sensitivity, and with which no separation in sensors for the determination of the position and in sensors for the determination of external fields is necessary.

While embodiments and applications of this invention have been shown and described, it would be apparent to those skilled in the art having the benefit of this disclosure that many more modifications than mentioned above are possible without departing from the inventive concepts herein. The invention, therefore, is not to be restricted except in the spirit of the appended claims and their equivalents. 

What is claimed is:
 1. Method for determining the position φ of a position indicator of a position measuring system, wherein the position measuring system comprises the position indicator and a number of at least r=2 sensors each delivering a sensor signal a_(j) with j=1 to r, the method having a calibration phase and a normal operation, wherein in normal operation the position φ of the position indicator is determined with the steps: A) detecting the sensor signals a₁ to a_(r) of the sensors and forming a signal vector {right arrow over (a)}=(a₁, . . . , a_(r)), B) forming a measurement vector {right arrow over (q)} according to a predetermined linear function g to {right arrow over (g)}=g({right arrow over (a)}), wherein the measurement vector {right arrow over (q)} has a number of n components, with 2≦n≦r, C) calculating a vector {right arrow over (p)}=M·{right arrow over (q)}, wherein {right arrow over (p)} is a 2-component vector with the components p₁ and p₂ and M is a 2×n matrix, and D) determining the position e by means of a predetermined function ƒ({right arrow over (p)}) to φ=ƒ({right arrow over (p)}), wherein the predetermined function ƒ({right arrow over (p)})=ƒ(p₁,p₂) is based on the equations p₁=cos(φ) and p₂=sin(φ), and wherein the matrix M is determined in a calibration phase either with the steps: detecting the sensor signals a_(k1) to a_(kr) of the sensors for a predetermined number of k=1 to m≧3 positions of the position indicator, forming m signal vectors {right arrow over (a)}k=(a_(k1), . . . , a_(kr)) from the sensor signals with k=1 to m, forming m measurement vectors {right arrow over (q)}_(k) with k=1 to m according to the predetermined function g to {right arrow over (q)}_(k)=g({right arrow over (a)}k), forming a matrix Q according to ${Q = \begin{pmatrix} {\overset{->}{q}}_{1} \\ \vdots \\ {\overset{->}{q}}_{m} \end{pmatrix}},$  and determining the matrix M from the matrix Q in such a way, that the transformation {right arrow over (p)}=M·{right arrow over (q)} is a linear mapping, which maps each of the m measurement vectors {right arrow over (q)}_(k) on a 2-component vector {right arrow over (p)}_(k), wherein the tips of the vectors {right arrow over (p)}_(k) substantially lie on a circle or on a circular arc, or by means of a bootstrap-method with the steps: setting an index i=0 and configuring a predetermined starting matrix M⁽⁰⁾, detecting and storing a plurality of measurement vectors {right arrow over (q)} according to the steps A and B of the normal operation, repeatedly executing the following steps: setting the matrix M to M=M^((i)) increasing the index i by 1, calculating the positions φ of the stored measurement vectors {right arrow over (q)} according to the steps C and D of the normal operation, selecting a number of k=1 to m of the stored measurement vectors {right arrow over (q)}_(k), the corresponding positions φ_(k) of which fulfill a predetermined criterion, forming a matrix Q according to ${Q = \begin{pmatrix} {\overset{->}{q}}_{1} \\ \vdots \\ {\overset{->}{q}}_{m} \end{pmatrix}},$ determining a new matrix M^((i)) from the matrix Q in such a way, that the transformation {right arrow over (p)}=M·{right arrow over (q)} is a linear mapping, which maps each of the m selected measurement vectors {right arrow over (q)}_(k) on a 2-component vector {right arrow over (p)}_(k), until the matrix M^((i)) within a predetermined tolerance no longer differs from the matrix M^((i−1)), whereby the tips of the vectors {right arrow over (p)}_(k) with increasing index i substantially come to lie on a circle or on a circular arc.
 2. Method according to claim 1, wherein the position measuring system comprises at least r=3 sensors.
 3. Method according to claim 1, wherein the determining of the matrix M comprises the following steps: determining two matrices U and B, which fulfill the equation Q=U·B, wherein U is a m×2 matrix with U^(T)·U=

, and B is a symmetrical, positive definite 2×2 matrix, by forming an auxiliary matrix D to D=Q^(T)·Q, calculating the matrix B according to B=D^(1/2), forming an auxiliary matrix H to H=B⁻¹, forming the matrix M to M=NH, wherein $N = \begin{pmatrix} d & 0 \\ 0 & 1 \end{pmatrix}$  is an auxiliary matrix, whose parameter d has the value 1 or −1.
 4. Method according to claim 2, wherein the determining of the matrix M comprises the following steps: determining two matrices U and B, which fulfill the equation Q=U·B, wherein U is a m×2 matrix with U^(T)·U=

, and B is a symmetrical, positive definite 2×2 matrix, by forming an auxiliary matrix D to D=Q^(T)·Q, calculating the matrix B according to B=D^(1/2), forming an auxiliary matrix H to H=B⁻¹, forming the matrix M to M=NH, wherein $N = \begin{pmatrix} d & 0 \\ 0 & 1 \end{pmatrix}$  is an auxiliary matrix, whose parameter d has the value 1 or −1.
 5. Method according to claim 1, wherein the determining of the matrix M comprises the following steps: determining three matrices U, S and V^(T), which fulfill the equation Q=U·S·V^(T), wherein U is an orthogonal m×m matrix, V^(T) is in orthogonal n×n matrix and S is a positive semidefinite m×n diagonal matrix, the diagonal elements of which are sorted according to decreasing size, forming sub-matrices Ũ, {tilde over (S)} and {tilde over (V)}^(T) in dependence of a predetermined parameter h, which fulfills the condition 2≦h≦l, wherein l≦min(m,n) designate the number of diagonal elements of S which are different from zero and wherein Ũ is that m×h sub-matrix of the matrix U, which is formed by discarding the columns h+1 to m, {tilde over (S)} is that h×h sub-matrix of the matrix S, which is formed by discarding the rows h+1 to m and discarding the columns h+1 to n, and {tilde over (V)}^(T) is that h×n sub-matrix of the matrix V^(T), which is formed by discarding the rows h+1 to n. if n≧2 forming an auxiliary matrix H from the elements of the sub-matrices Ũ, {tilde over (S)} and {tilde over (V)}^(T), or in the special case of n=2 alternatively forming the auxiliary matrix H to H={tilde over (S)}⁻¹{tilde over (V)}^(T), and forming the matrix M zu M=NH, wherein $N = \begin{pmatrix} d & 0 \\ 0 & 1 \end{pmatrix}$  is an auxiliary matrix, whose parameter d has the value 1 or −1.
 6. Method according to claim 5, wherein the forming of the auxiliary matrix H from elements of the sub-matrices Ũ, {tilde over (S)} and {tilde over (V)}^(T) comprises the following steps: forming an auxiliary matrix W, if the measurement vectors are present in an arbitrary order with the steps: calculating values α_(i)=a tan 2(U₁₂U_(i1)−U₁₁U_(i2), U₁₂U_(i2)+U₁₁U_(i1)) for i=1 to n, calculating the permutation k(i), which sorts the α_(i)-values according to their size, so that α_(k(i))>α_(k(i−1)) for i=2 . . . , m, and forming the elements of the auxiliary matrix W according to $W_{i,{k{(j)}}} = \left\{ {\begin{matrix} {{\sin \left( \frac{2{\pi \left( {j - 1} \right)}}{m} \right)},{i = 1}} \\ {{\cos \left( \frac{2{\pi \left( {j - 1} \right)}}{m} \right)},{i = 2}} \end{matrix},} \right.$ for j=1 to m, if the m positions of the position indicator are distributed over approximately the entire range from 0 to 2π, or $W_{i,{k{(j)}}} = \left\{ {\begin{matrix} {{\sin \left( \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m} \right)},{i = 1}} \\ {{\cos \left( \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m} \right)},{i = 2}} \end{matrix},} \right.$ for j=1 to m if the m positions of the position indicator occupy only a circular arc with φ∈[φ_(min),φ_(max)], or $W_{i,{k{(j)}}} = \left\{ {\begin{matrix} {{\sin \left( {e \cdot \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m}} \right)},{i = 1}} \\ {{\cos \left( {e \cdot \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m}} \right)},{i = 2}} \end{matrix},} \right.$ for j=1 to m, if the m positions of the position indicator occupy only a circular arc with φ∈[φ_(min),φ_(max)] and the position indicator has a multipole character, wherein the parameter e characterizes the multipole character, or if the measurement vectors {right arrow over (q)}_(k) with k=1, . . . , m are sorted in ascending or descending order according to their corresponding position φ_(k), forming the elements of the auxiliary matrix W with the auxiliary quantity ƒ=1 for ascending and ƒ=−1 for descending order according to $W_{ij} = \left\{ {\begin{matrix} {{\sin \left( {f \cdot \frac{2{\pi \left( {j - 1} \right)}}{m}} \right)},{i = 1}} \\ {{\cos \left( {f \cdot \frac{2{\pi \left( {j - 1} \right)}}{m}} \right)},{i = 2}} \end{matrix},} \right.$ for j=1 to m if the m positions of the position indicator are distributed over approximately the entire range from 0 to 2π, or $W_{ij} = \left\{ {\begin{matrix} {{\sin \left( {f \cdot \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m}} \right)},{i = 1}} \\ {{\cos \left( {f \cdot \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m}} \right)},{i = 2}} \end{matrix},} \right.$ for j=1 to m if the m positions of the position indicator occupy only a circular arc with φ∈[φ_(min),φ_(max)], or $W_{ij} = \left\{ {\begin{matrix} {{\sin \left( {{fe} \cdot \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m}} \right)},{i = 1}} \\ {{\cos \left( {{fe} \cdot \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m}} \right)},{i = 2}} \end{matrix},} \right.$ for j=1 to n if the m positions of the position indicator occupy only a circular arc with φ∈[φ_(min),φ_(max)] and the position indicator has a multipole character, wherein the parameter e characterizes the multipole character, forming an auxiliary matrix T to T=WŨ, forming the auxiliary matrix H to H=T{tilde over (S)}⁻¹{tilde over (V)}^(T).
 7. Method according to claim 2, wherein the determining of the matrix M comprises the following steps: determining three matrices U, S and V^(T), which fulfill the equation Q=U·S·V^(T), wherein U is an orthogonal m×m matrix, V^(T) is in orthogonal n×n matrix and S is a positive semidefinite m×n diagonal matrix, the diagonal elements of which are sorted according to decreasing size, forming sub-matrices Ũ, {tilde over (S)} and {tilde over (V)}^(T) in dependence of a predetermined parameter h, which fulfills the condition 2≦h≦l, wherein l≦min(m,n) designate the number of diagonal elements of S which are different from zero and wherein Ũ is that m×h sub-matrix of the matrix U, which is formed by discarding the columns h+1 to m, {tilde over (S)} is that h×h sub-matrix of the matrix S, which is formed by discarding the rows h+1 to m and discarding the columns h+1 to n, and {tilde over (V)}^(T) is that h×n sub-matrix of the matrix V^(T), which is formed by discarding the rows h+1 to n. if n≧2 forming an auxiliary matrix H from the elements of the sub-matrices Ũ, {tilde over (S)} and {tilde over (V)}^(T), or in the special case of n=2 alternatively forming the auxiliary matrix H to H={tilde over (S)}⁻¹{tilde over (V)}^(T), and forming the matrix M zu M=NH, wherein $N = \begin{pmatrix} d & 0 \\ 0 & 1 \end{pmatrix}$  is an auxiliary matrix, whose parameter d has the value 1 or −1.
 8. Method according to claim 7, wherein the forming of the auxiliary matrix H from elements of the sub-matrices Ũ, {tilde over (S)} and {tilde over (V)}^(T) comprises the following steps: forming an auxiliary matrix W, if the measurement vectors are present in an arbitrary order with the steps: calculating values α_(i)=a tan 2(U₁₂, U_(i1)−U₁₁U_(i2), U₁₂U_(i2)+U₁₁U_(i1)) for i=1 to m, calculating the permutation k(i), which sorts the α_(i)-values according to their size, so that α_(k(i))>α_(k(i−1)) for i=2, . . . , m, and forming the elements of the auxiliary matrix W according to $W_{i,{k{(j)}}} = \left\{ {\begin{matrix} {{\sin \left( \frac{2{\pi \left( {j - 1} \right)}}{m} \right)},{i = 1}} \\ {{\cos \left( \frac{2{\pi \left( {j - 1} \right)}}{m} \right)},{i = 2}} \end{matrix},} \right.$ for j=1 to m, if the m positions of the position indicator are distributed over approximately the entire range from 0 to 2π, or $W_{i,{k{(j)}}} = \left\{ {\begin{matrix} {{\sin \left( \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m} \right)},{i = 1}} \\ {{\cos \left( \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m} \right)},{i = 2}} \end{matrix},} \right.$ for j=1 to m if the m positions of the position indicator occupy only a circular arc with φ∈[φ_(min),φ_(max)], or $W_{i,{k{(j)}}} = \left\{ {\begin{matrix} {{\sin \left( {e \cdot \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m}} \right)},{i = 1}} \\ {{\cos \left( {e \cdot \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m}} \right)},{i = 2}} \end{matrix},} \right.$ for j=1 to m, if the m positions of the position indicator occupy only a circular arc with φ∈[φ_(min), φ_(max)] and the position indicator has a multipole character, wherein the parameter e characterizes the multipole character, or if the measurement vectors {right arrow over (q)}_(k) with k=1, . . . , m are sorted in ascending or descending order according to their corresponding position φ_(k), forming the elements of the auxiliary matrix W with the auxiliary quantity ƒ=1 for ascending and ƒ=−1 for descending order according to $W_{ij} = \left\{ {\begin{matrix} {{\sin \left( {f \cdot \frac{2{\pi \left( {j - 1} \right)}}{m}} \right)},{i = 1}} \\ {{\cos \left( {f \cdot \frac{2{\pi \left( {j - 1} \right)}}{m}} \right)},{i = 2}} \end{matrix},} \right.$ for j=1 to m, if the m positions of the position indicator are distributed over approximately the entire range from 0 to 2π, or $W_{ij} = \left\{ {\begin{matrix} {{\sin \left( {f \cdot \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m}} \right)},{i = 1}} \\ {{\cos \left( {f \cdot \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m}} \right)},{i = 2}} \end{matrix},} \right.$ for j=1 to m if the m positions of the position indicator occupy only a circular arc with φ∈[φ_(min),φ_(max)], or $W_{ij} = \left\{ {\begin{matrix} {{\sin \left( {{fe} \cdot \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m}} \right)},{i = 1}} \\ {{\cos \left( {{fe} \cdot \frac{\left( {\varphi_{\max} - \varphi_{\min}} \right)\left( {j - 1} \right)}{m}} \right)},{i = 2}} \end{matrix},} \right.$ for j=1 to m if the m positions of the position indicator occupy only a circular arc with φ∈[φ_(min),φ_(max)] and the position indicator has a multipole character, wherein the parameter e characterizes the multipole character, forming an auxiliary matrix T to T=WŨ, forming the auxiliary matrix H to H=T{tilde over (S)}⁻¹{tilde over (V)}^(T).
 9. Method according to claim 2, wherein ƒ(p₁,p₂)=a tan 2(p₂,p₁)−φ₀, wherein φ₀ is a quantity which characterizes a zero position of the position measuring system.
 10. Method according to claim 7, wherein ƒ(p₁,p₂)=a tan 2(p₂,p₁)−φ₀, wherein φ₀ is a quantity which characterizes a zero position of the position measuring system.
 11. Method according to claim 8, wherein ƒ(p₁,p₂)=a tan 2(p₂,p₁)−φ₀, wherein φ₀ is a quantity which characterizes a zero position of the position measuring system.
 12. Method according to claim 2, wherein the function g({right arrow over (a)}) is one of the following functions: 1) g({right arrow over (a)})={right arrow over (a)}, 2) g({right arrow over (a)})={right arrow over (a)}−{right arrow over (a)}₀, 3) g({right arrow over (a)})=P·{right arrow over (a)}, or 4) g({right arrow over (a)})=P·({right arrow over (a)}−{right arrow over (a)}₀), wherein {right arrow over (a)}₀ designates an offset vector and P a projection matrix.
 13. Method according to claim 7, wherein the function g(a) is one of the following functions: 1) g({right arrow over (a)})={right arrow over (a)}, 2) g({right arrow over (a)})={right arrow over (a)}−{right arrow over (a)}₀, 3) g({right arrow over (a)})=P·{right arrow over (a)}, or 4) g({right arrow over (a)})=P·({right arrow over (a)}−{right arrow over (a)}₀), wherein {right arrow over (a)}₀ designates an offset vector and P a projection matrix.
 14. Method according to claim 13, wherein g({right arrow over (a)})=P·{right arrow over (a)} or g({right arrow over (a)})=P·({right arrow over (a)}−{right arrow over (a)}₀) and wherein the projection matrix P is determined in the calibration phase, either in absence of the position indicator with the steps: executing the following steps for a predetermined number of different external fields v=1 to b≧1, wherein the number b fulfills the condition b≦r−2: applying the external field, detecting the sensor signals a₁ to a_(r) and forming a signal vector {right arrow over (a)}=(a₁, . . . , a_(r)), and forming an external field vector {right arrow over (ƒ)}_(v)={right arrow over (a)} or, given the case, of an offset compensated external field vector {right arrow over (ƒ)}_(v)={right arrow over (a)}−{right arrow over (a)}₀, or in presence of the position indicator with the steps: moving the position indicator in a fixed position, detecting the sensor signals a₁ to a_(r) and storing as vector {right arrow over (a)}_(F)=(a₁, . . . , a_(r)), executing the following steps for a predetermined number of different external fields v=1 to b≧1, wherein the number b fulfills the condition b≦r−2: applying the external field, detecting the sensor signals a₁ to a_(r) and forming a signal vector {right arrow over (a)}=(a₁, . . . , a_(r)), and forming an external field vector {right arrow over (ƒ)}_(v)={right arrow over (a)}−{right arrow over (a)}_(F) or, given the case, of an offset compensated external field vector {right arrow over (ƒ)}_(v)={right arrow over (a)}−{right arrow over (a)}_(F)−{right arrow over (a)}₀, and both in the presence and absence of the position indicator either with the further steps: forming of c orthonormal vectors {right arrow over (g)}₁, . . . , {right arrow over (g)}_(c) of the vector space spanned by the external field vectors {right arrow over (ƒ)}₁, . . . , {right arrow over (ƒ)}_(b), wherein c≦b, complementing the c vectors {right arrow over (g)}₁, . . . , {right arrow over (g)}_(c) with n=r−c further vectors {right arrow over (g)}_(c+1), . . . , {right arrow over (g)}_(c+n) to a complete, orthonormal basis of the r-dimensional vector space of the signal vectors, and forming the projection matrix P according to $P = {\begin{pmatrix} {\overset{->}{g}}_{c + 1} \\ \vdots \\ {\overset{->}{g}}_{c + n} \end{pmatrix}.}$ or with the further steps: forming an orthonormal basis of c vectors {right arrow over (g)}₁, . . . , {right arrow over (g)}_(c) of the vector space spanned by the external field vectors {right arrow over (ƒ)}₁, . . . , {right arrow over (ƒ)}_(b), and forming the projection matrix P according to ${P = { - {\sum\limits_{i = 1}^{c}{{\overset{->}{g}}_{i} \otimes {\overset{->}{g}}_{i}}}}},$  wherein {right arrow over (g)}_(i)

{right arrow over (g)}_(i) are the tensor products of the vectors {right arrow over (g)}_(i).
 15. Method according to claim 8, wherein the function g(d) is one of the following functions: 1) g({right arrow over (a)})={right arrow over (a)}, 2) g({right arrow over (a)})={right arrow over (a)}−{right arrow over (a)}₀, 3) g({right arrow over (a)})=P·{right arrow over (a)}, or 4) g({right arrow over (a)})=P·({right arrow over (a)}−{right arrow over (a)}₀), wherein {right arrow over (a)}₀ designates an offset vector and P a projection matrix.
 16. Method according to claim 15, wherein g({right arrow over (a)})=P·{right arrow over (a)} or g({right arrow over (a)})=P·({right arrow over (a)}−{right arrow over (a)}₀) and wherein the projection matrix P is determined in the calibration phase, either in absence of the position indicator with the steps: executing the following steps for a predetermined number of different external fields v=1 to b≧1, wherein the number b fulfills the condition b≦r−2: applying the external field, detecting the sensor signals a₁ to a_(r) and forming a signal vector {right arrow over (a)}=(a₁, . . . , a_(r)), and forming an external field vector {right arrow over (ƒ)}_(v)={right arrow over (a)} or, given the case, of an offset compensated external field vector {right arrow over (ƒ)}_(v)={right arrow over (a)}−{right arrow over (a)}₀, or in presence of the position indicator with the steps: moving the position indicator in a fixed position, detecting the sensor signals a₁ to a_(r) and storing as vector {right arrow over (a)}_(F)=(a₁, . . . , a_(r)), executing the following steps for a predetermined number of different external fields v=1 to b≧1, wherein the number b fulfills the condition b≦r−2: applying the external field, detecting the sensor signals a₁ to a_(r) and forming a signal vector {right arrow over (a)}=(a₁, . . . , a_(r)), and forming an external field vector {right arrow over (ƒ)}_(v)={right arrow over (a)}−{right arrow over (a)}_(F) or, given the case, of an offset compensated external field vector {right arrow over (ƒ)}_(v)={right arrow over (a)}−{right arrow over (a)}_(F)−{right arrow over (a)}₀, and both in the presence and absence of the position indicator either with the further steps: forming of c orthonormal vectors {right arrow over (g)}₁, . . . , {right arrow over (g)}_(c) of the vector space spanned by the external field vectors {right arrow over (ƒ)}₁, . . . , {right arrow over (ƒ)}_(b), wherein c≦b, complementing the c vectors {right arrow over (g)}₁, . . . , {right arrow over (g)}_(c) with n=r−c further vectors {right arrow over (g)}_(c+1), . . . , {right arrow over (g)}_(c+n) to a complete, orthonormal basis of the r-dimensional vector space of the signal vectors, and forming the projection matrix P according to $P = {\begin{pmatrix} {\overset{->}{g}}_{c + 1} \\ \vdots \\ {\overset{->}{g}}_{c + n} \end{pmatrix}.}$ or with the further steps: forming an orthonormal basis of c vectors {right arrow over (g)}₁, . . . , {right arrow over (g)}_(c) of the vector space spanned by the external field vectors {right arrow over (ƒ)}₁, . . . , {right arrow over (ƒ)}_(b), and forming the projection matrix P according to ${P = { - {\sum\limits_{i = 1}^{c}{{\overset{->}{g}}_{i} \otimes {\overset{->}{g}}_{i}}}}},$  wherein {right arrow over (g)}_(i)

{right arrow over (g)}_(i) are the tensor products of the vectors {right arrow over (g)}_(i).
 17. Method according to claim 5, wherein a matrix C is formed in the calibration phase as C=

−{tilde over (V)}·{tilde over (V)} ^(T), and wherein in normal operation a check value v is formed with the steps forming a vector {right arrow over (d)}=C{right arrow over (q)}, calculating the check value v from the vector {right arrow over (d)}, and triggering a predetermined action, if the check value v exceeds a predetermined value.
 18. Method according to claim 6, wherein a matrix C is formed in the calibration phase as C=

−{tilde over (V)}·{tilde over (V)} ^(T), and wherein in normal operation a check value v is formed with the steps forming a vector {right arrow over (d)}=C{right arrow over (q)}, calculating the check value v from the vector {right arrow over (d)}, and triggering a predetermined action, if the check value v exceeds a predetermined value.
 19. Method according to claim 1, further comprising in running operation capturing and buffering a plurality of measurement vectors {right arrow over (q)} and of the corresponding positions φ which have been calculated according to the steps C and D of the normal operation, and performing a recalibration by means of the steps: selecting a number of buffered measurement vectors {right arrow over (q)}, the corresponding positions φ of which fulfill a predetermined criterion, rebuilding the matrix Q from the selected measurement vectors {right arrow over (q)}, and rebuilding the matrix M.
 20. Method for determining the position φ of a position indicator of a position measuring system, wherein the position measuring system comprises the position indicator and a number of r sensors each delivering a sensor signal a_(j) with j=1 to r, the method having a calibration phase and a normal operation, wherein in normal operation the position φ of the position indicator is determined with the steps A) detecting the sensor signals a₁ to a_(r) of the sensors and forming a signal vector {right arrow over (a)}=(a₁, . . . , a_(r)), B) forming a measurement vector {right arrow over (q)} according to a predetermined linear function g to {right arrow over (q)}=g({right arrow over (a)}), wherein the measurement vector {right arrow over (q)} has a number of n=2 components, wherein the function g({right arrow over (a)}) is one of the following functions: 1) g({right arrow over (a)})=P·{right arrow over (a)}, or 2) g({right arrow over (a)})=P·({right arrow over (a)}−{right arrow over (a)}₀), wherein {right arrow over (a)}₀ designates an offset vector and P a projection matrix, and C) determining the position φ by means of a predetermined function ƒ({right arrow over (q)}) to φ=ƒ({right arrow over (q)}), and wherein in the calibration phase the projection matrix P is determined with the steps: either in absence of the position indicator with the steps: executing the following steps for a predetermined number of different external fields v=1 to b=r−2: applying the external field, detecting the sensor signals a₁ to a_(r) and forming a signal vector {right arrow over (a)}=(a₁, . . . , a_(r)), and forming an external field vector {right arrow over (ƒ)}_(v)={right arrow over (a)} or, given the case, of an offset compensated external field vector {right arrow over (ƒ)}_(v)={right arrow over (a)}−{right arrow over (a)}₀, or in presence of the position indicator with the steps: moving the position indicator in a fixed position, detecting the sensor signals a₁ to a_(r) and storing as vector {right arrow over (a)}_(F)=(a₁, . . . , a_(r)), executing the following steps for a predetermined number of different external fields v=1 to b≧1, wherein the number b fulfills the condition b=r−2: applying the external field, detecting the sensor signals a₁ to a_(r) of the sensors and forming a signal vector {right arrow over (a)}=(a₁, . . . , a_(r)), and forming an external field vector {right arrow over (ƒ)}_(v)={right arrow over (a)}−{right arrow over (a)}_(F) or, given the case, of an offset compensated external field vector {right arrow over (ƒ)}_(v)={right arrow over (a)}−{right arrow over (a)}_(F)−{right arrow over (a)}₀, and both in the presence and absence of the position indicator with the further steps: forming of c orthonormal vectors {right arrow over (g)}₁, . . . , {right arrow over (g)}_(c) of the vector space spanned by the external field vectors {right arrow over (ƒ)}₁, . . . , {right arrow over (ƒ)}_(b), wherein c=b, complementing the c vectors {right arrow over (g)}₁, . . . , {right arrow over (g)}_(c) with two further vectors {right arrow over (g)}_(c+1), {right arrow over (g)}_(c+2) to a complete, orthonormal basis of the r-dimensional vector space of the signal vectors, and forming the projection matrix P according to $P = {\begin{pmatrix} {\overset{->}{g}}_{c + 1} \\ {\overset{->}{g}}_{c + 2} \end{pmatrix}.}$ 